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Abstract 

The growth and viriaUzation of spherical top-hat fluctuations, in coupled Dark Energy mod- 
els, causes segregation between Dark Matter (DM) and baryons, as the gravitational infall into 
the potential well proceeds more slowly for the baryons than for DM. As a consequence, after 
attaining their turn-around and before full virialization, halos have outer layers rich of baryons. 
Accordingly, a natural ambiguity exists on the definition of the virial density contrast. In fact, 
when the outer baryon layers infall onto the DM-richer core, they carry with them DM materials 
outside the original fluctuation; hence, no time exists when all materials originally belonging to 
the fluctuation - and only them- have virialized. Baryon-DM segregation can have various as- 
trophysical consequences on different length-scales. The smallest halos may loose up to 50 % of 
the original baryonic contents and become hardly visible. Subhalos in cluster-size halos may loose 
much baryonic materials, which could then be observed as intra-cluster light. Isolated halos, in 
general, can be expected to have a baryon component richer than the cosmological proportions, 
due to the cosmic enrichement of baryons lost in small halo encounters. 

PACS numbers: 98.80.-k, 98.65.-r 
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INTRODUCTION 



Dark Energy (DE) was first required by SNIa data, but a flat cosmology with Qm — 0.27, 
h ~ 0.7 and ^ 0.045 was soon favored by CMB and LSS observations {Qm,b'- matter, 
baryon density parameters; h: Hubble parameter in units of 100 km/s/Mpc; CMB: cosmic 
microwave background radiation; LSS: large scale structure). The nature of Dark Energy is 
one of the main puzzles of cosmology. If it arises from false vacuum, its pressure/energy- 
density ratio w = Pde/pde is -1. For vacuum to yield DE, however, a severe fine tuning at 
the end of the electroweak transition is required. Otherwise, DE could be a scalar field (f) 
self-interacting through a potential V{(f)) Q], so yielding 

PDE = Pk,DE + Pp,DE = 9/2a^ + V{(j)), PdE = Pk,DE - Pp,DE (l) 

provided that pk/V <^ 1/2; then —1/3 ^ w > —1. Here, the background metric reads 
ds"^ = a'^l—dr'^ + dxidx^) {i = 1,..,3) and dots yield differentiation with respect to the 
conformal time r. This kind of DE is dubbed dynamical DE (dDE) or quintessence and a 
large deal of work has been done on it (see, e.g., and references therein). 

At variance from Dark Matter (DM), DE fluctuations tend to fade on scales signifi- 
cantly smaller than the horizon. Furthermore, while non-gravitational interactions between 
baryons and DE can be safely excluded by experimental data, constraints on DM-DE inter- 
action - as on most DM properties - can only be derived from cosmological observations. 
In turn, such interaction could ease the cosmic coincidence problem j4|, i.e. that DM and 
DE densities, after being different by orders of magnitude for most of the cosmic history, 
approach equal values in today's Universe. In a number of papers, constraints on coupling, 
coming from CBR and LSS data, were discussed 0,0,0]. 

In this note we focus on coupled DE (cDE) models and study the growth of a spheri- 
cal top-hat fluctuation. This is done assuming a SUGRA ^ or a RP |2] self-interaction 
potential: 

SUGRA V{(f)) = (A"+70°) exp(47r0V"^J) (2) 

RP V{(P) = (A°+V0") (3) 

(mp = G^^/^: Planck mass) as self-interaction potential. Once Ude assigned, either a or 
the energy scale A can be freely chosen. 



The background continuity equations for DE and DM then read 

+ 2(d/a)0 + aV,^ = ^levrG/S/^aVc , pc + 3(d/a)pc = IQtiG /?,l3pc^ , (4) 

while radiation and baryons obey standard equations. The interaction between DE and DM 
is parametrized by (3 which, all through this paper, is assumed constant and, in particular, 
not to depend on (j). Using WMAP data on T anisotropies, a limit (3 <~ 0.3 has been 
established. A more stringent limit on /?, found for the potential by studying halo 
profiles in cDE simulations p], does not apply to this case. 

In this paper we explored different values of /3, taking A = 10^ GeV and fi^ = 0.30, 
= 0.05 (model I); where significant, this was compared with a cosmology with VLm = 0.25, 
Qf, = 0.04 (model II); h = 0.7 was taken for both models. Full detail will be given for model 

I with SUGRA potential. 

This model provides a fair fit to CMB data, even performing slightly better than ACDM 
0]. Model II and RP potential are then considered just to explore how results depend on 
the choice of cosmological parameters and potential. 

Accordingly, results for Model I and SUGRA will be compared with Model II and SUGRA 
and/or Model I and RP. Some quantitative outputs will be however provided also for Model 

II and RP. 

The evolution of a spherical top-hat fiuctuation, in the linear and non-linear regimes, can 
be found analytically in a standard CDM model. Here the density contrast at virialization 
Aj, ~ 180. Although assuming spherical symmetry, which is not the geometry of growing 
fiuctuations, this treatment, associated with a Press & Schechter ^9] - or similar Q| - 
formulation, predicts observational and numerical mass functions with a high degree of 
approximation. The value of A^, is also used to extract virialized systems from galaxy 
samples and n-body simulations. 



The analysis of the spherically simmetric growth has been extended to ACDM mode 
ll| and other models with uncoupled DE, where it was used with similar success 
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In these cosmologies, the simple system of equations yielding spherical evolution requires a 
numerical solution. 

Models with DM-DE coupling are also considered by Q| but they do not take into 
account the different dynamics between coupled and uncoupled components as well as their 
bias (see eq.ljHl) below). 
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Infact, in the case of cDE, the dynamics of baryons and DM are different. An effective 
Newtonian theory can be safely used to describe their physical evolution, when radiative 
materials are dynamically irrelevant and for scales well below the horizon scale, both for 
small and non-linear fluctuations, as was done, e.g., by j^. 

In the newtonian limit, the DM-DE interaction causes: (i) DM particle masses to vary 
oc exp(— C0); (ii) the gravitational constant between DM particles to become G* = 7G. 
Here 

C = ^le^rG/S (3 , 7 = 1 + 4/3V3 . (5) 

Taking this into account, at the initial time Tj, we assume that both DM and baryons have 
a top-hat fluctuation of identical radius RTH,ii expanding together with the Hubble flow. 
The linear theory prescribes that fluctuation amplitudes in baryons and DM start from 
a ratio 

J, ~ 37 l]c + 4/3X/i ■ ^ ' 

Here VLc is the DM density parameter, /x = {6c,b/Sc,b)/{o,/a), X = ^JAti/S ^ / {rripa / a) . Due 
to the different interaction strength, as soon as a non-linear regime is reached, the size i?* 
of the baryon perturbation begins to exceed the size i?^ of the DM perturbation. Hence, a 
part of the baryons initially within Rth,! leak out from R'^, so that DM does not feel the 
gravity of all baryons, while baryons above R'^ feel the gravity also of initially unperturbed 
DM layers. As a consequence, above R'^, the baryon component deviates from a top-hat 
geometry, while a secondary perturbation in DM arises also above i?^ itself. 

After reaching maximum expansion, contraction will also start at different times for 
different components and layers. Then, inner layers will approach virialization before outer 
layers, whose later fall-out shall however perturb their virial equilibrium. This already 
outlines that the onset of virial equilibrium is a complex process. 

Furthermore, when the external baryon layers fall-out onto the virialized core, richer of 
DM, they are accompanied by DM materials originally outside the top-hat, perturbed by 
baryon over-expansion. 

The time when the greatest amount of materials, originally belonging to the fluctuation, 
are in virial equilibrium occurs when the DM top-hat has virialized, together with the 
baryon fraction still below R'^. A large deal of baryonic materials are then still falling out. 
But, when they will accrete onto the DM-richer core, they will not be alone, carrying with 
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them originally alien materials. There will be no discontinuity in the fall-out process when 
all original baryons are back. The infall of outer materials just attains then a steady rate. 

In order to follow the dynamical evolution of a systems where each layer feels a substance- 
dependent force, a set of concentric shells, granting a sufficient radial resolution, needs to 
be considered. In the next section we shall provide the equations of motion for Lagrangian 
shell radii and discuss the whole expansion and recontraction dynamics. More details on 
the virialization process are provided in Section 3. Section 4 has a technical content and 
is devoted to working out the Newtonian regime from General Relativistic equations. In 
Section 5 a brief discussion on the relevance of baryon-DM segregation, on various length- 
scales, is given. Final conclusions are drawn in Section 6. 

TIME EVOLUTION OF CONCENTRIC SHELLS 

In order to write dynamical equations, we shall make use of the comoving radii 6„ = R\^/a 
and Cn = Rn/(^ for the n-th baryon or DM shell, respectively. The masses of layers depend 
on the Rn distribution at the initial time. Let M^''^ be the mass of the n-th layer for each 
component and let M^''^ be the mass for the same layer and component in the absence of 
the perturbation. Accordingly, we can define the mass excess ^M^''^ = M^''^ — M^''^, for the 
ra-th layer. While the masses of baryonic shells keep constant while they expand and then 
recontract, according to eq. (jH) it is 

M;:(r) = M;:(r,„)exp{-C[0(r) - 0(r,„)]} • (7) 

Similarly, let AM'''%< r) = M^'^(< r)~Mb.c{< r) be the mass excess within the physical 
radius R = ar and let then be 

AM^ = AM\< hn) + AM\< hn) , AM^ = AM\< c„) + 7AM'^(< c„) , (8) 

taking into account both the mass variation (eq. Hj) and that G* ^ G. Here AM^'^(< r) 
are masses, approximately given by the sum of all the 5M^^'^ for which 6„, < r. AM^''^, 
instead, are gravity sources, given by masses corrected - somewhere- by suitable factors, 
different for DM and baryons. 

Problems due to the systematic discrepancy between the comoving radii 6„ and c„ can 
be cured by increasing the resolution and taking into account only the fraction of the last 
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shell below the radius considered. The shell dynamics can then be described through the 
equations 



eqs. (jH) and the Friedman equation. 

The eqs. (jSj) are analogous for DM and baryons. Until 6„ ~ c„, it is AM^ < AM^ (eq. EI), 
for any (3 > 0. Hence, the gravitational push felt by DM layers is stronger. The extra 



derivative. As a consequence, the DM comoving radius c„ decreases more rapidly than the 
baryon comoving radius 6„ (or increases more slowly than R^) and the ra-th baryon 
radius may gradually exceed the [n + l)-th DM shell, etc. . As a consequence, the sign of 
AM^ — AM^ can even invert and one can evaluate for which value of — 0i this occurs. 
Once 0(r) is known, also the time when this occurs can be found. Until then, however, DM 
fluctuations expand more slowly and mostly reach their turn-around point earlier, while 
baryons gradually leak out from the fluctuation bulk. 

This is visible in the Figure ^ for model I. Here the time dependence of a sample of 
radii R!^'^ is shown. The greatest radii shown are the top-hat radii. Solid and dotted lines 
yield the a dependence of the ratios R^/Ri and R^/Ri respectively, for (3 = 0, 0.1, 0.2, 0.3 . 
These plots are obtained from the solution of eqs. Q down to the collapse of the last baryon 
shell, assumed to occur at a = 1 (today). An extreme case (/3 = 0.7) is also shown in Fig. |21 
The actual number of radii R^'^ used depends on the level of precision wanted and radii 
do not need to be equi-spaced. The requirement fulfilled here is that, doubling the radius 
number in the critical regions, where profiles are expected to suffer a stronger deformation, 
causes a shift in all outputs never exceeding 1 %. To this aim, up to ~ 1000 radii had to 
be used. If more precision is wanted, it can be easily achieved at the expenses of using a 
greater computer time. With a single processor 1.5 Gflops PC, the numerical program takes 
~ 90 minutes to run a single case, with ~ 1000 radii. 

Full collapse is not expected in the physical CclSG, clS the Rn decrease shall stop when 
virialization is attained and, in the next section, the virialization condition will be discussed. 
Figures El IH El however, assume that the spherical growth stops when all DM originally in 
the top-hat, and the baryons kept inside it, virialize. 

In Figure El the effects of the spherical growth are shown in the limit of no coupling 





which are derived in Sec. 4. These equations are to be integrated together with the first of 



term C0c„ adds to this push. In fact, the comoving variable Cn has however a negative 
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FIG. 1: Evolution of a sample of baryon and DM layer radii, extending up to the top-hat radius, 
in model I, for couplings /? up to 0.3 . 

{P = 0) for model I. This plots allows to outline an effect that is usually disregarded, i.e., 
that also materials outside from the top-hat are perturbed by its growth. Their density, 
in fact, increases in respect to average. As shown in Figure ^ the turn-around occurs for 
a ~ 0.6. At such time, the density enhancement of exterior layers extends up to a radius 
approximately double of the top-hat radius. When the virialization radius is reached by 
DM and inner baryons (a ~ 0.92), the perturbed sphere is even wider and, at twice the 
top-hat radius, the density enhancement is ~ 2.4. These effects, however, do not bear a 
great relevance in the (3 = case, where layer sharing is unessential, as growth is self-similar 
over all shells inside the top-hat. 

In Figure EJ the case (3 = 0.1 is considered. Already at a = 0.2, well before the turn- 
around, the slope of the top-hat boundary, for baryons, is no longer vertical. At a = 0.6 
(approximately turn-around), not only the baryon boundary is bent, but a similar effect is 
visible also for DM, where the density increase at -R > Rth, is however greater than for 
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FIG. 2: Evolution of a sample of baryon and DM layer radii in model I, for /? = 0.7 . Although 
such coupling is somehow large, the plot allows the reader to see the detail of layer evolution. 

P = 0. The effect is even more pronounced at amr — 0.92. Figure El refers to /? = 0.3, where 
all above effects are present and stronger. 

The increased density of DM shells outside the top-hat is relevant, here, because it 
fastens the recoUapse of outer baryons. In turn, the enhanced baryon density, outside from 
DM top-hat, modifies the dynamics of external DM layers, as well. 

Let us also outline that the initial fluctuation amplitude is different in Figures ^ and 01 
EJ In the former case, it is tuned so that the outest layer fully recollapses at a = 1. In 
the latter case, instead, only the DM fluctuation would have fully recoUapsed at a = 1. 

Another point which ought to be noticed in the same Figures is that the only density 
discontinuity occurs at the boundary of the DM top-hat and is due to the DM component 
only. Baryon density begin to decline there, but is not dicontinuous. The baryon excess, 
due to the initial top-hat, makes them intially denser than DM a,t R > The density 
ratio inverts at some outer radius and the baryon density decrease is faster than for outer 
DM layers. The radius i?^, where the most external baryon top-hat shell lays, is visible in 
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FIG. 3: Density profiles at different o values for uncoupled DE model. Solid (dotted) lines refer 
to DM (baryons), In the absence of coupling the two profiles have the same shape. The shell 
treatment, in this case, in unnecessary; notice, however, the density increase outside the top-hat. 

most plots, as a discontinuity of the density slope. It however lays even below 2 R'^, in some 
profiles. 

As outlined before, the most apparent physical effect is the outflow of baryon layers from 
the DM top-hat. The outflown baryon fraction increases with a. In the top panel of Figure 
such fraction is shown for (3 varying from 0.1 to 0.3 in the case of SUGRA potential. 
For a ~ 0.92, i.e. when DM and inner baryons have attained their virialization radius, the 
fraction of baryons which have leaked out from the fluctuation is ~ 20 % for (3 = 0.1 and 
reaches ~ 58 % for (3 = 0.3 . These values increase by an additional 10 % in RP case. This 
is shown in the top panel of Figure [7| The lower panel of both Figs. IHl and [7| shows that the 
effect depends mildly on the scale A and increases when A is greater. 
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FIG. 4: Density profiles at different a values for (3 = 0.1 model. Notice the progressive deformation 
of the baryon profile (dotted lines) in respect of the DM profile (solid line). 
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FIG. 5: Density profiles at different a values for /3 = 0.3 model. The distorsion of the baryon 
profile (dotted line), in this case, is already noticeable at turn-around and, around virialization, 
the baryon profile has almost lost its initial top-hat shape. 
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FIG. 6: Baryon fraction outside top-hat approaching viriahzation (SUGRA model). 




FIG. 7: Baryon fraction outside top-hat approaching viriahzation (RP model). 
VIRIALIZATION 

As already outlined, in Figures |5 El the DM profile preserves its initial top-hat config- 
uration which, on the contrary, is gradually distorted in the baryon component. In order 
to set the viriahzation time, therefore, one must first state to which shells it refers. An 
obvious option amounts to relating viriahzation to top-hat materials which are not polluted 
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by accretion of external layers. This definition (option I), however, leaves out a large deal 
of baryons, which will accrete later on, although accompanied by alien DM. Another choice 
(option II) is then to wait until all baryons are allowed to virialize, regardless of the intruding 
DM, which alterates the original proportions. Any intermediate choice is also legitimate. 
For instance, one could refer to the time when the total mass reaching its virialization radius 
equates the total mass intially inside the top-hat, regardless of the DM/baryon ratio, which 
will be however finally increased. 

Most of this discussion does however risk to be non-physical. The study of the evolution 
of a spherical top-hat is significant when applied to a realistic cosmic environment. We 
shall take this point of view in Section 5. Here below we shall define that virialization has 
occurred when DM and inner baryons would have fully recoUapsed (option I). 

During fluctuation evolution, the potential energies both of DM and baryons will be made 
of three terms, arising by self-interaction, mutual interaction, and interaction with the DE 
field: 

U^R) = V'iR) + U''\R) + f/"'^^(i?) = 47r ^ dr r"^ Pc{r) \^<^{r) + ^^(r) + ^ DE{r)] (10) 

Jo 

U\R) = U''\R) + U'^^R) + [/''•^^(i?) = An ("^ dr r'' p^ir) [^^(r) + <^,{r) + ^oEir)] (11) 

Jo 

Here the sum on the layers is replaced by a volume integral, extended to a sphere of radius 
R. While 

Att Att 
"^bir) = -— Gp,(r)r2 , vl>^^(r) = -—GpnEiry , (12) 

in both expressions, a subtle difference exists for "if^- In U'' we have simply 

Att 

^c(r) = --Gp^iry , (13) 

as for the other components, but this expression is different in If^, where it reads 

An 

^c(r) = --{lG[p, - Pc(r)] + Gp^y , (14) 

Pc being the background DM density. The different dynamical effect of background and DM 
fluctuation arises from the different ways how its interaction with DE is treated. Energy 
exchanges between DM and DE, for the background, are accounted for by the r.h.s. terms 
in eqs. In this case no newtonian approximation was possible and was made. For DM 
fluctuation, instead, the effects of DM-DE exchanges are described by a correction to the 
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gravitational constant G, becoming G* = 7G, which adds to the dependence of DM density 
on (f). This ought to be taken into account in the fluctuation evolution, as is done in eq. (|T^. 
Together with the potential energies, the kinetic energies 

T'^iR) = \j^ dVp, f2 , T\R) = 1 1^ dVp, (15) 

for DM and baryons within i?, due to the collapse motion, ought to be considered. The 
virialization condition then reads 

2T{R) = RdU{R)/dR (16) 

Here T{R) = T''{R)+T\R) and U{R) = U^{R) + U\R). If i? is selected within the DM top- 
hat, Pb,c,DE do not depend on r (such dependence would become important if virialization 
were defined taking into account layers outside the DM top-hat) and then integrals (fTUj). 
(fTT|) and (fT^ can be performed obtaining 

5 it 5 

5 it 5 

T^(i?) = 1 - M'R^ (19) 
2 5 

where the relation r/r = R/R is used to calculate T'^{R). Note that this relation is not 
valid for T^{R) because different baryon layers have different growth rates. Kinetic energy 
for baryons is then obtained by 

T\R) = Y.T'n = Y.\M'SRif (20) 

Here the sum is extended on all i?^ < R. Notice that, because of DE-DM interaction yielding 
DM particle mass variation, we cannot require energy conservation for DM and baryons 
between turn-around and virialization. However, energy conservation in the spherical growth 
is violated, also in cosmologies with uncoupled dynamical DE, once we assume full DE 
homogeneity. This problem as well as the possibility to relax the homogeneity assumption 
for DE are discussed by 

In Figure IHl we report the /^-dependence of the density contrast at virialization, both for 
models I and II, both for SUGRA and RP potential. 
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FIG. 8: Density contrast at virialization in cDE models. 
THE NEWTONIAN REGIME 

This section debates the passage from the General Relativistic treatment to the Newto- 
nian Regime equations, for cDE theories. The discussion partially resumes arguments made 
by 15 1 and Q], whose outcomes are used to work out eqs. ((Tj). The reader, who is not 
interested in deriving these equations, can skip this section. 

When fluctuations are considered, in the Newtonian gauge the metric reads 



ds'^ = a2(r)[-(l + 2ip)dT^ + (1 - 2ip)dxidx% (21) 

ip being the gravitational potential. No anisotropic stresses are considered here. As usual, 
fluctuations are expanded in components of wavenumber k and let be \ = l-L/k [Ti. = d/a). 
Dealing with fluctuations well after recombination, radiation is neglected. Baryon (DM) 
density fluctuations are dubbed 6b (5c); DE field fluctuations 6(f) and (DM or baryon) velocity 
fields \c,b are described by the variables 



^ = y^'50, ^ ^^^^ 
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Without restriction on the DE potential V{(f)) one can define the functions 



/ = rV3/167rG ln(F/K), h=^^ + f^ = 0^ + 2/ + A ; (23) 

{Vo is a reference value of the potential); it is also useful to define X"^ = 47rG'0^/37i^ (as in 
the introduction) and Y'^ = 8TTGV{(j)) a^/37Y^. The gravitational and (p field equations then 
read 

i/j = -^X\nbSb + ^cSc + GX^ + 2X'p' -2Y^fi^) , ij' = 3x^ - ij , (24) 

h' 

f" + (2 + ttt)^' + - 12Xy. + 4^X + 2Y\hv - M) = P^dSc + 2^), (25) 
rt 

keeping just the lowest order terms in A, to obtain the Newtonian limit; if DE kinetic (and/or 
potential) energy substantially contributes to the expansion source, X (and/or Y) is 0{1). 

In the Newtonian limit, we must also neglect the derivatives of ip, averaging out the 
oscillations of ip and the potential term f2Y'^ip, requiring that A << (j^Y)^^ (remind that 
Y is 0{1)). Furthermore, in eq. (j2Sl), the metric potential ip (oc A^) can also be neglected. 
Accordingly, eq. ()24|) and (^^j) become 

^ = -^A2(n,56 + nA), , A-V-/9^^c4. (26) 

(the former one is the Poisson equation). Then, from the stress-energy pseudo-conservation 
T^^ = 0, one can derive a couple of equations telling us how 6cfi and 6*^6 depend on a: 

5," = -5,'{1 + ^ - 2/3X) + ^(1 + ^/3')nA + l^bSb, (27) 

5b" = -5b' (1 + %) + l{^c5c + ^b5b). (28) 
rt I 

= -e^l + 2(5X) - \{\ + ^/32)1]A - %^b5b. , (29) 
rt Z 6 Z 

Ob ' = -Obil + ^) - l{^c5c + ^b5b) ; (30) 
rt Z 

here ' yields differentiation with respect to a = In a. Assuming Vtb « Vt^. and putting 
5c oc e-f ^^"^^^ and 5b = b5c with b = cost from the eqs.(j2Zl), (EH), we obtain the bias factor 
(jni). The acceleration of a single DM or baryon particle of mass m^fe can be instead derived 
from eqs. (jSHI); (jSDI)- Let us set it in the void, at a distance r from the origin, where a DM (or 
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baryon) particle of mass Mc (or Mb) is set, and let us remind that, while the usual scaling 
Pb oc holds, it is 

= p,,a-3e-^(<^-<^°), pM. = M,,a-3e-^('^-<^o)(5(0), (31) 

because of the DE-DM coupling; here the subscript o indicates values at the present time t^, 
(it is tto = 1). We can then assign to each DM particle a varying mass Mc{(f)) = Moce~'"^'^~'^°^ 
Then, owing to eq. (jHH), and assuming that the density of the particle widely exceeds the 
background density, it is 

= '-^^ = ^MM)m , ^.Sb = = ^MbKO), (32) 

Per iK'^a Per Sn'^a 

[pcr is the critical density and 6 is the Dirac distribution). Reminding that V ■ Vc^^ = O'^'^'Ti. 
and using the ordinary (not conformal) time, eq. (j27p yields 

V ■ V, = -H{1 - 2(3X) V ■ Ve - 47rGa"2 {lM,{(j)) + Mb) 6{0) (33) 

(dots yield differentiation in respect to ordinary time and H = a /a). Taking into account 
that the acceleration is radial, as the attracting particles lie at the origin, it will be 



J d^rV ■ V = 477 J dr d{r'^v)/dr = Arcr'^v. 



Accordingly, the radial acceleration of a DM particle read 

..^-,l-2,A-)i/v..n-^^-^f, (34) 

(n is a unit vector in the radial direction; R = ar). 

Repeating the calculation for a baryon we get immediately the result 

Vb = -H^b ■ n — — (35) 

In the presence of a full spherical symmetry, Vb,c ■ n = Vb^c and using b = Vb and c = r^, it is 
easy to derive eqs. (0). 

For the sake of completeness, here below, we write also the equations for the physical 
DM and baryons radii which easily follow from eqs. ()34|) and ()35|) . While for baryons we 
have the usual Friedman-like equations 

Ait 

Ri = - [pe + Pb + Pde{1 + 'iw)] Ri , (36) 

for DM we have 

Rl = C^K - C^HR^^ - [p, + 7(Pc - Pc) +Pb + Pde{1 + Sw)] R'^. (37) 
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DISCUSSION 



In this work we studied the evolution of a top-hat fluctuation, assuming it spherical and 
that its growth proceeds without external perturbations. Neither assumption is realistic. 
In spite of that, it is known that the results of a spherical unperturbed growth, obtained 
considering only one spatial coordinate, yield a significant insight into the evolution of 
fluctuations in the 3-dimensional space. 

It seems clear, e.g., that baryon outflow from original density fluctuations is not peculiar 
for spherical overdensities. The fraction of baryons laying outside, at the various growth 
stages, on the contrary, can only be considered as an estimate. 

If external perturbations are completely negligible, the general outcome is that the final 
virialized system is however richer of DM. This is true indipendently of the way how virial- 
ization is defined. If we define it referring to DM and inner baryons only, the contribution of 
the outer baryon layers is clearly missing. But, even if we wait and define virialization when 
all baryons, initially belongong to the fluctuations, have accreted, a substantial amount of 
DM, previously external to the top-hat, has also accreted by then. In both cases the final 
DM/baryon ratio has apparently increased and this is true also if we stop at any intermediate 
stage. 

If we look closer at the quantitative details, it appears that the DM excess is greatest 
when we refer to DM and inner baryons only and becomes gradually smaller when we include 
further layers. The actual DM/baryon ratio of real systems will therefore depend on the 
stage up to which we take accretion into account. In turn, this is linked to the second 
assumption, that the growth occurs without external perturbations. Steady accretion can 
be assumed to end when perturbations can no longer be neglected. Let us return on this 
point below. 

It would however be fallacious to deduce that the DM/baryon ratio, in all stabilized 
systems, exceeds cosmological proportions. When a DM richer system is formed, baryons 
are not destroyed, but enrich the cosmic medium. Systems forming later will therefore start 
their growth in a baryon richer environment. When their formation proceeds unperturbed 
until the last baryon layer has accreted, the residual DM excess can be smaller than the 
initial baryon excess and the system can be even richer of baryons than cosmic proportions 
prescribe. 
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The above arguments seem therefore to indicate that a gradual decrease of the average 
DM/baryon ratio will occur when going from smaller to greater systems. Large system form 
later and suffer less important perturbations. In cDE theories, therefore, galaxy clusters can 
be expected to preserve or exceed the cosmological baryon content. 

Inside large halos, yielding galaxy clusters, smaller subhalos continue to grow. Inside a 
cluster, the interaction of systems is frequent and outer layers are likely to be stripped. We 
can therefore expect, in cDE theories, that galaxies in clusters are DM richer, while a large 
deal of baryon materials are stripped by tidal forces in close encounters, presumably form 
baryon rich filaments and finally enrich the intracluster medium (ICM). 

Observations show that a hot intracluster plasma with temperatures T ~ 10''-10^K emits 
over a broad band from EUV to X-ray, from the potential well of galaxy clusters. Most 
observed radiation features are accounted assuming that the emission process is thermal 
)reemstrahlung, plus lines associated with the metals in ICM. Many observations in EUV 
lyj l and soft X-ray however claim an excess emission. Large n-body simulations (see, 
e.g., have been used to try to locate the origin of this excess, which apparently arise 
from filaments at T ~ 10^-lO^K . Simulations apparently succeed in explaining the EUV 
and soft X-ray excess, but fail to locate it in filaments. In turn, filaments could be a natural 
product of outer baryon stripping. 

It is clearly premature to suggest that baryon-DM segregation helps to reach an agree- 
ment between data and simulation outputs. There are however good grounds to assess that 
the whole picture of ICM emission is modified by segregation. 

The existence of a large amount of baryons in ICM is also testified by optical data 
based, e.g., on planetary nebulae observations (see, e.g^2^ and references therein). Again, 
simulations meet data with some reserves (see, e.g., |2l| and references therein). Clearly, 
baryon stripping from forming sub-halos would substantially modify model predictions. 
In particular, stars in ICM could turn out to be not so older than the average galactic 
populations. 

When the smallest galactic systems are then considered, a clear predictions of cDE theo- 
ries is that they can be significantly baryon poor. Another prediction is that their individual 
baryon contents will be variable, depending on the frequency and timing of close encounters 
with greater structures. It seems however rather clear that a galaxy satellite, before being 
dynamically settled, can easily meet bigger objects and suffer tidal stripping of outer layers. 
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FIG. 9: Upper panel: Dependence on (3 of the virial mass. Lower panel: Dependence on [3 of the 
virialization radius. 



221 and 



In the recent years there has been a wide discussion on observed galaxy satellites. 
[23I 1 showed that, in a ACDM cosmology, the number of predicted satellite galaxies exceeds 
the observed number in the local group by more than one order of magnitude. Predictions for 



ACDM are confirmed also for dynamical DE j2^. Various mechanisms were then proposed, 
by a number of authors, to reduce the abundance of substructure j^] or to suppress star 



26| 



formation in small clumps via astrophysical mechanisms, making them dark 

Because substructure can modify the fluxes of lensed images relative to those predicted by 
smooth lens models, gravitational lensing can be used to detect and constrain populations 
of DM c 
authors 



umps in galactic halos j27|. Different methods are, however, proposed by several 

We cannot state here that the baryon depletion of very small halos, envisaged here, is the 
solution of the above problem. The question however deserves to be more deeply examined 
through suitable techniques. 
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CONCLUSIONS 



In this work we studied the spherical growth of top-hat fluctuations in cDE theories. 
Our findings concern: (i) The ambiguity of the definition of halo virialization, within such 
context which, in turn, could cause some difficulty in comparing simulations outputs or data 
with Press & Schechter or similar predictions, (ii) The baryon-DM segregation during the 
spherical growth. 

The (i) effect is critical if one wants to use results on the gravitational growth of a 
spherical fluctuation to perform predictions on the mass function for large bound systems, 
like galaxy clusters. In fact, the mass enclosed in the virialized region, indipendently of its 
definition, does not coincide with the mass enclosed in a fluctuation evolving according to 
the linear Jeans' equations. 

In upper panel of Figure IHl we give the ratio between the virialized mass for various j3 
values and the virialized mass for (3 = 0, coinciding with the mass in a linearly evolving 
fluctuation. The shift arises both from a different virial density contrast (see Fig. |HI) and a 
different virial radius (Figure El lower panel). 

The construction of a cluster mass function should take into account these outputs. In a 
forthcoming paper, we shall deal with this question in full detail. 

Various consequences of the (ii) effect, on different scales, were also envisaged. On very 
large scales it leads to predict some baryon enrichments of large clusters. On intermediate 
scales it does interfere with a number of tentative matches between simulations and obser- 
vations. On small scales it could offer a way out to the empasse arising from the scarsity of 
observed galaxy satellites in the Local Group. 

Quite in general, this analysis shows that cDE theories can open interesting perspectives 
in the treatment of a number of cosmological problems. In turn, observations can be expected 
to put precise limits on the strength of the possible DM-DE coupling. 
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